Class D spectral peak in Majorana quantum wires 
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Proximity coupled spin-orbit quantum wires purportedly support midgap Majorana states at 
critical points. We show that in the presence of disorder these systems generate a second bandcenter 
anomaly, which is of different physical origin but shares key characteristics with the Majorana state: 
it is narrow in width, insensitive to magnetic fields, carries unit spectral weight, and is rigidly tied to 
the band center. Depending on the parity of the number of subgap quasiparticle states, a Majorana 
mode does or does not coexist with the impurity peak. The strong 'entanglement' between the two 
phenomena may hinder an unambiguous detection of the Majorana by spectroscopic techniques. 



More than seven decades ago, the Majorana fermion 
was postulated as a neutral variant of the Dirac parti- 
cle pp. While the Majorana has never been observed 
in the canonical environments of particle physics, recent 
proposals [2H1] suggest a realization as a boundary state 
of one dimensional topological superconductors. The per- 
spective to realize, and potentially apply topologically 
protected Majorana particles in condensed matter set- 
tings, has sparked a wave of theoretical [5HT5] and ex- 
perimental [TSUI!?] activity. 

Recently, three experimental groups have reported sig- 
natures of Majorana states in semiconductor quantum 
wires. In these experiments, proximity coupled semicon- 
ductor quantum wire purportedly hosting a Majorana 
particle are probed by tunneling spectroscopy. Evidence 
for the presence of the particle is drawn from the obser- 
vation of a zero energy spectral peak. It has been verified 
that the presence of the peak hinges on parametric condi- 
tions necessary for the formation of a Majorana particle. 

In this Letter we argue that the quasi one-dimensional 
quantum wires presently under consideration are prone 
to the formation of a second zero energy peak, caused by 
strong midgap quantum interference. The latter struc- 
ture, here called (class D) spectral peak for brevity [20] . 
shares key universal signatures with the Majorana peak, 
to the extent that disentangling the two phenomena may 
be difficult: it is (i) tied to zero energy, (ii) affected by 
temperature, but not (iii) by magnetic field, (iv) relies 
on parametric conditions similar to those required by 
the Majorana, (v) carries unit spectral weight (a sin- 
gle quasi 'state'), (vi) is independent of other midgap 
structures, such as Kondo resonances or Andreev bound 
states, shows (vii) sensitivity to the parity of channel 
number, however, unlike the Majorana, (viii) relies on 
the presence of moderately weak disorder. 

Qualitative discussion and main results. — To start 
with, let us qualitatively discuss the physics underlying 
the spectral peak. The spin orbit and proximity coupled 
quantum wires presently under discussion form a variant 
of a spinless superconductor, more precisely a system of 
class D [3T]. At the critical configuration separating a 
topologically trivial and nontrivial state (which will be 
the habitat of the Majorana fermion if parameters are 




FIG. 1. a) Phase coherent scattering process p — > h — > p, 
b) low energy diffusion mode formed by succession of such 
processes, c) cartoon of a higher order quantum interference 
process involving self-interacting diffusion modes. 

tuned in a space-like manner) the superconductor goes 
gapless. The presence of a Majorana particle signifies 
in a band center anomaly of the quasiparticle density of 
states (DoS). Quantum mechanically, the DoS at energy 
e (relative to the system's chemical potential, //) is given 
by p(e) — — —J dx Im tr(G + (x, x; e)), where the Green 
function G + (x,y,e) = (x\(e + iO + p,o- ph — H)^ 1 ^), ag h is 
a Pauli matrix in the superconductor particle hole (ph) 
space, and 'tr' is a trace over that space. Semiclassically, 
the matrix elements G pp (x, x; e) (or G hh (x, x; e)) may be 
interpreted as return amplitudes of quasi-particles, in 
an environment governed by the Bogoliubov-de Genncs 
Hamiltonian H of the system. Much like in a gapless 
superconductor |22j . the DoS in the proximity quantum 
wire is affected by the scattering of low energy excitations 
off spatial fluctuations of an order parameter amplitude, 
A. To understand the physics of this mechanism |21) . 
imagine a quasi-particle of energy e emanating from the 
point x. The particle has the option to scatter off the 
order parameter A into a quasi-hole of energy — e, which 
by a second scattering process may get converted back 
to a particle returning to the point of origin. For low ex- 
citation energies, the particle and hole stretches involved 
in this process interfere constructively to give a robust 
(scattering phase insensitive) contribution to the return 
amplitude. Summation over repeated ph conversion pro- 
cesses generates an effective diffusion mode [ST] similar to 
the diffuson and Cooperon modes of conventional disor- 
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dered metals. The increase of the effective quasi-particle 
gap away from the critical point confines the support of 
this mode to a small region in space - a disordered gap- 
less superconductor 'quantum dot' - where it scales as 
~ 8/ (27re) = s _1 , 8 being the dot's effective level spacing. 
In disordered electronic systems, singularities of this type 
are generally cut off by the nonlinear 'self-interaction' of 
diffusion modes (weak localization turning into strong 
localization, or level repulsion being prominent examples 
of such nonlinearities). In the present context, similar 
nonlinearities lead us to the result 

p( s ) sin(s) „ e 

iVeven: t^- = l+ — ±-L, s = 2tt- 
Po s 8 

Nodd: M = 1 -^. + S (JL), (1) 

where N is the number of channels in the system, 5 the 
average quasiparticle level spacing, po = 8 , and the 8- 
function is the contribution of the Majorana state. Mul- 
tiple quantum interference multiplies the single diffusion 
mode by a phase factor ±sin(s), where the overall sign 
depends on the channel number parity. The ensuing pro- 
files of the band center DoS are shown in the inset to 
Fig. [2] Notice the superficially similar structures for even 
and odd channel numbers, while only in the latter case 
a genuine Majorana mode is present. For even channel 
numbers, a narrow width (~ 8) peak is solely formed by 
class D diffusion modes. In either case, the anomalous 
contribution to the DoS integrates to the spectral weight 
J de(p(e) — po) = |. This means that tunneling spec- 
troscopy limited in resolution to scales larger than 8 may 
not be able to unambiguously identify the Majorana state 
in situations where quantum interference is effective. An 
external loss of coherence, e.g. due to quasiparticle inter- 
actions, will suppress the class D peak, although a power 
law correction ~ | e | 1 to the density of states p(e) above 
an effective cutoff energy should remain visible. Finally 
notice how the phenomenon hinges on the presence of 
disorder: impurities render the low energy dynamics dif- 
fusive, and they provide the basis for the formation of 
an 'impurity band' of quasiparticles centered around the 
gap closing point in energy and space. However, we will 
argue below that even weak scattering t _1 > 8 suffices 
to generate the effect. The rest of the paper is devoted 
to a derivation of Eq. ([!]). We will also discuss physical 
bounds on the spatial extension of the impurity quantum 
dot, the relevant disorder strengths, and other system pa- 
rameters. 

Model. — Majorana states can be realized in different 
variants of class D topological superconductors. For defi- 
niteness, we here consider the case of a proximity coupled 
helical liquid [llj subject to a smoothly varying mag- 
netic field. However, we will argue below that our main 
conclusions are not rigidly tied to this setup. A helical 
liquid is formed by a system of left- and right-moving 
fermions, ipLrf = "0L an d ipR,l = ipR carrying spin up 
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FIG. 2. Schematic profile of spectral density. The profile of 
the band center anomaly depends on the parity of the number 
of channels, N. Only for odd N the system supports a genuine 
Majorana level which (here represented as a broadened 5- 
singularity.) 

and down, respectively. (In the semiconductor setting, 
these modes form at the intersection of fermion bands 
shifted in momentum by strong spin-orbit interaction.) 
The coupling to a wire axis Zeeman field, B, and a prox- 
imity s-wave order parameter then generates the effective 
Hamiltonian H Q = J2a=i I dx (H a K + H% +H a A ), where 
K a K = Ec=l,r r c is cVFd x r c , Kb = BV L V R + h-c, 
H% — Atplipft + h.c, vf is the Fermi energy, and 
sc = (+/—)! for C — L/R. We have also introduced 
an index a — 1,...,N, which accounts for the option 
of multiple bands tp^ crossing the Fermi energy. Differ- 
ences in the Fermi velocities of the bands are absorbed 
in a rescaling of the fields. Ignoring the effects of band 
coupling by spin orbit interaction, the effects of static 
disorder and interaction, etc. Hq has the status of a 
null theory describing the formation of a critical regime 
at B ~ A: introducing new state vectors ('Majorana 

basis' 1 ! n° = 1 ( -^H~^« ^ n 1 = 1 ( «/>R-"?R ^ a „A 

basis ), V _ y^[^ L _^ L ),fl - ^ [ ) , and 

assuming reality of the order parameter, A ~ B, the 
Hamiltonian Hq assumes the form 




where a summation over channel indices is implicit and 
the Pauli matrices act in Li?-space. Notice that the ma- 
trix structure sandwiched by the 77's is antisymmetric, 
which is the defining condition of a class D (no symme- 
tries other than particle-hole symmetry) superconductor. 
The above representation makes the formation of a low- 
(77 ) and a high-energy sector (r/ 1 ) manifest. In the ab- 
sence of perturbations, rj°' a represent Majorana modes 
which go massless at the critical point A e ff = A — B = 0. 
Throughout, we will assume a hierarchy of energy scales 
M = B + A ^> A c ff, and on this basis we now turn 
to the discussion of perturbations, 8H, to H . Scat- 
tering between the low- and the high-energy sector - 
generated by fluctuations of the chemical potential, or 
imaginary contributions to the order parameter - affects 
the low energy theory by contributions of 0(SH 2 /M)) 
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and will be neglected throughout. By contrast, real 
fluctuations of the order parameter, or channel depen- 
dent differences in the field strength (indirectly gener- 
ated by inter-channel spin-orbit coupling) couple within 
the low energy sector and have to be kept. We model 
such fluctuations in terms of a phenomenological ansatz 
A e ff = Ao(ir) + W + V, where the first term describes 
a ramping of system parameters through a critical point 
at x = 0, where A (0) = (cf. Fig. ||, W = {W ab } is 
a matrix containing deterministic contributions to the 
inter-channel coupling, and V = {V ab (x)} represents 
static disorder, here assumed to be Gaussian distributed, 
(V ab (x)V a ' b '(x')) = S aa 'd bb 'S(x - x')7 2 . Throughout we 
will limit our discussion to channels for which the low en- 
ergy physics around x = is disorder dominated, i.e. for 
which the random potential V > W masks the determin- 
istic contribution. (If this condition cannot be satisfied 
by a multiplet of channels, we are down to N — I .[23 ) 

a-model. — Our focus throughout will be on regions 
close to x = where the gap amplitude |A(x)| < V 
is small enough for a class D impurity quantum dot 
to form. Such structures universally exhibit a zero en- 
ergy spectral peak which we aim to explore quantita- 
tively. To this end, we consider the supersymmetric func- 
tional integral J D(£, £)exp(i J dx£(e + — H)£), where 
H = —ivpd x a3 — A c ff(x)<72 is the low energy (/i = 0) 
Hamiltonian, and £ = {£c,a} an d £ = {£c,a} are inte- 
gration variables with complex (a = b) and Grassmann 
valued (a = f) components. The DoS can be obtained 
from the above functional by differentiation w.r.t. suit- 
ably introduced sources, which we make not explicit for 
notational simplicity. 




FIG. 3. Top: a sweep of an effective gap parameter generates 
a critical point at which the clean system goes gapless (cf. 
schematic dispersion relations shown in the insets.) Bottom: 
in the presence of disorder an impurity band of semicircular 
local average DoS forms. 

We next average the generating functional over V, 
decouple the ensuing term quartic in £ by Hubbard- 
Stratonovich transformation, and integrate over £, 
to obtain the action S[X,X] = ^2 / dxstr(XX) + 



^strln(e + T 3 + Ag + A^c^ — Ho), where 'str' is the super- 
trace [H] and H = H I 



v=o' 



Here we have introduced a 



new 2-component 'cc' space whose sectors ~ (tp, y ) con- 
tains states, ip, and their Majorana-basis charge conju- 



gates, ip T , resp. The Pauli-matrix, T3 acting in that space 
reflects the different transformation behavior of H and 
e + , (e + —H) T = —(—£ + —H). The four component super- 
matrices X = {X ta t i a >}, and X = XVj f are Hubbard- 
Stratonovich fields, where the Pauli matrix in bf space 
ensures convergence, t, if are indices in cc-space, and 
Xs,A — \{X±X). Following standard strategies [23] . we 
process the above action by a variational approach. Em- 
ploying the polar ansatz, X = X = iXTAT^ 1 , where T 
are rotation matrices to be discussed momentarily and A 
is a diagonal matrix carrying real matrix elements ±1, we 
find that the modulus A = X(x) > is fixed by the vari- 



ational equation X(x) = 7 2 Im(a;|(iA 
upon substitution of Hq becomes 



H ) 



E 
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dp 

2^ (v FP ) 2 + A a2 + A 2 



E 



VA a2 + A 2 



which, 



= 1. 



Here we have assumed the matrix A c ff| y o = {A a (5 } 
diagonalized, and varying slowly enough in space for the 
functions X(x) to adiabatically follow. Inspection of the 
equation shows that solutions exist if |A a | < Nj 2 /vp on 
average over a. In the limiting case of negligible channel 
dependence, A a ~ A , we readily obtain 



X(x) = 




1/2 



-Al(x) 



where r _1 = 2Nj 2 /vp is the golden rule scattering rate 
of the unperturbed spin-orbit quantum wire. The inter- 
pretation of Gq (e) = (e + — Ho + iX)^ 1 as an averaged 
Green function then yields v = 2tvqX for the (position 
dependent) DoS per unit length at saddle point level, 
where vq = N/(irvp). From this quantity, we gener- 
ate the average DoS as po = §,-. >Q dxv(x), where the 
integral is over the region of support of the local DoS, 
|Ao(a;)| < l/2r (cf. Fig. [3]) In passing we note that for 
large energies |e| 3> t _1 , the density of states assumes a 
form reminiscent of a smeared BCS profile to be discussed 
in more detail elsewhere [2"5] . 

Turning to the saddle point matrix A, we note that, as 
with the standard cr-model |24j . the eigenvalue structure 
of the bosonic sector Abb = T3 is fixed by the pole signa- 
ture of the Green function. However, in the fermionic sec- 
tor we have a choice, Ag = ±T3, an ambiguity intimately 
related to the formation of a Majorana peak in the 
limit of low energies, e — > 0, spatially uniform rotations 
A — > TAT -1 = Q leave the action invariant. The symme- 
tries of our class D theory require [27] Tbb € Sp(2)/U(l), 
while T ff e 0(2)/U(l) ~ Z 2 . This latter identity means 
that in the fermionic sector we have only two discrete 
symmetry transformations, Tg = 1, and Tg = n, where 
the latter acts by an exchange Tgr^Tg 1 = —T3 of the 
diagonal saddle points. There are no continuous trans- 
formations T interpolating between Ai and A2, which 
means that the saddle point manifold, M. = Mi U A^2 
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splits into two disjoint components, where Ai s , s — 1, 2 is 
generated by action of all transformations T with Tff = 1 
on the saddle point A = Ai = t$ and A = A2 = T3 ® 03, 
resp. 

We finally substitute slowly fluctuating configurations 
X — > iXQ into the action and expand in small symmetry 
breaking contributions ~ e + Q and dQ. The expansion 
follows standard procedures [22j [24] and generates the 
effective action 

S[Q] = J dx 7n ^tr(D(x)dQdQ + 4eT 3 Q) + S T>aI [Q} ) 

where S par [Q] = |strln(A — Hq) is the Pfaffian of the 
model Hamiltonian at saddle point level, and D — 
^ S a (a 2 +a° 2 ) 3 / 2 Pl & y s tne T0 ^ e °f an effective diffusion 
coefficient, which in the isotropic limit A a ~ Ao simpli- 
fies to D ~ (vpr) 2 X. Notice how in the disorder domi- 
nates regions |A| <C A, D ~ v\t reduces to the diffusion 
constant of a metallic system, while D — > at the bound- 
aries A — > of the metallic domain. The structure of the 
action above indicates that fluctuations of Q generate the 
diffusive soft modes discussed in connection with Fig. [l] 
Quasiparticle 'quantum dot'. - For energies e < 
(D) /L lower than the minimum excitation energy of 
spatially varying Q-fluctuations (L is the extension of the 
metallic domain A > |Ao|, and (D) is the typical value 
of the diffusion coefficient), the action is dominated by 
homogeneous modes Q = const. In this regime, the sys- 
tem truly behaves like a 'quantum dot', and the action 
collapses to 

S[Q] = ^-tr(Qr 3 ) + 5 par [Q]. (2) 

Here, S par is a topological term which in the present 
context assumes the form S pax [Q] = S paT [Ax] = and 
S pax [Q] = S par [A 2 ] = trln(G+(0)-Go (0)) for Q e Mi, 2 
resp. Using the particle-hole symmetry of the spectrum, 
the latter expression evaluates to S^ar^] = — i^Vtot 
where AT to t = J^a^a, ana - ^ a * s the total number of 
states carried by mode a. While this number does not 
carry physical meaning by itself, exp(— S pai ) = (— ) Nt ° t 
responds only to the parity of N to t- Now, robust argu- 
ments show that N a = (2fc+l) is generally odd 28 . This 
means that (— ) JVtot = (— ) N is susceptible to the parity 
of the number of channels. The DoS resulting from Eq. 

computed in Ref. [H] by non-perturbative integra- 
tion over .Mi ,2, is given by Eq. ([I]). 

While Eq. ([IJ has been derived for a specific model, 
the high degree of universality of the result — depen- 
dence on no parameters other than the ratio e/S — testi- 
fies to wider applicability. Generally speaking, we expect 
the result to hold provided 5 < t _1 -C A, i.e. the ran- 
domness hybridizing a number of levels of the clean sys- 
tem around the critical point x ~ 0. This expectation is 
corroborated by the observation of similar spectral pro- 
files for phenomenological random matrix models |21j . 



(A straightforward numerical calculation indeed shows 
that peaks similar to (JlJ form already for two randomly 
coupled levels.) Finally, random spectra are notorious 
for their 'level rigidity', i.e. the weakness of sample-to- 
sample fluctuations. Signatures of the anomaly in the 
ensemble averaged global DoS are therefore expected to 
show in both the sample specific global and local DoS. 

Summarizing, we have shown that in the presence of 
even weak disorder, a Majorana peak in class D quan- 
tum wires will coexist with a spectral anomaly. The two 
phenomena share key universal features, are of compa- 
rable magnitude, and integrate to a spectral weight in- 
sensitive to the parity of the channel number. All this 
makes us speculate that disorder may be detrimental to 
an unambiguous observation of the Majorana particle by 
spectroscopic means. 

We thank F. von Oppen for discussions. Work sup- 
ported by SFB/TR12 of the Deutsche Forschungsgemein- 
schaft. 
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